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ABSTRACT 

Galactic black hole transients show many interesting phenomena during outburst decays. We present 
simultaneous X-ray {RXTE, Swift, and INTEGRAL), and optical/near- infrared (0/NIR) observations 
(SMARTS) of the X-ray transient XTE J1752-223 during its outburst decay in 2010. The multiwave- 
length observations over 150 days in 2010 cover the transition from soft to hard spectral state. We 
discuss the evolution of radio emission is with respect to the 0/NIR light curve which shows several 
flares. One of those flares is bright and long, starting about 60 days after the transition in X-ray 
timing properties. During this flare, the radio spectral index becomes harder. Other smaller flares 
occur along with the X-ray timing transition, and also right after the detection of the radio core. 
We discuss the significances of these flares. Furthermore, using the simultaneous broadband X-ray 
spectra including INTEGRAL, we find that a high energy cut-off with a folding energy near 250 keV 
is necessary around the time that the compact jet is forming. The broad band spectrum can be fitted 
equally well with a Comptonization model. In addition, using photoelectric absorption edges in the 
XMM-Newton RGS X-ray spectra and the extinction of red clump giants in the direction of the source, 
we find a lower limit on the distance of > 5 kpc. 

Keywords: black hole physics - X-rays:stars - accretion, accretion disks - binaries: close 



1. INTRODUCTION 

Galactic black hole transients (hereafter GBHTs) show 
distinct spectral and temporal properties during the 
whole outburst across all wavelengths. According to their 
X-ray spectral and timing properties, these systems are 
found mainly in the soft state (SS) or in the hard state 
(HS). Here, we provide a brief summary of properties 
of GBHTs in these states (see McClintock & Remillard 
2006; Belloni 2010, for details). In the SS, the X-ray 
emission is dominated by soft photons from a geomet- 
rically thin, optically thick disk, which is modeled by 
a multi-temperature blackbody (Makishima et al. 1986). 
The variability is low in the SS, often below the detection 
limit of RXTE for short observations. In the HS, the flux 
from the thin disk is weak, and the emission is dominated 
by a hard component, phenomenologically modeled with 
a power-law in the X-ray spectrum. The variability is 
strong (typically >20% rms amplitude), and the power 
spectra seen in this state can exhibit a variety of features 
such as quasi-periodic oscillations. There also exist in- 
termediate states (IS), usually observed when the source 
makes a transition between SS and HS. 

The emission properties in radio, optical and infrared 
bands also strongly depend on the X-ray spectal states. 
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The radio emission is quenched in the SS (Fender et al. 
1999; Corbel & Fender 2002; Russell et al. 2011). In 
the HS, radio observations indicate compact steady jets 
(Fender 2001), and synchrotron radiation from the com- 
pact jet has also been detected in optical, and near 
infrared (Buxton & Bailyn 2004; Kalemci et al. 2004; 
Migliari et al. 2007; Coriat et al. 2009; Buxton et al. 
2012; Kalemci et al. 2013). For XTE J1550-564, it has 
been reported that the direct jet synchrotron emission 
may dominate the entire X-ray emission in the hard state 
at 10-4 < LEdd < 10-3 (Russell et al. 2010). Simi- 
lar analysis indicated that the jet producing significant 
X-ray fiux at low luminosities is also a possibility for 
XTE J1752-223 (Russell et al. 2012). 

During outburst decays, these sources display state 
transitions from the SS to IS to HS. Since the jet emis- 
sion is quenched in the SS, multiwavelength observations 
of outburst decays allow us to investigate the X-ray spec- 
tral and timing properties required for jet formation and 
to constrain the jet's contribution to the X-ray emission. 

A possible effect of the jet on the X-ray spectral 
properties is the injection of non-thermal electrons in 
the corona, thereby hardening the X-ray spectrum. A 
cut-off in high energy spectrum of these sources is of- 
ten observed in the hard state, usually interpreted as 
a sign of thermal Comptonization. Sometimes the cut- 
off disappears in high energies after jets are observed 
as in 4U 1543-47 (Kalemci et al. 2005), and GRO 
J1655— 40 (Kalemci et al. 2006b; Caballero Garcia et al. 
2007)). However, there are sources that do not re- 
quire a cut-off for the entire outburst (H1743— 322 
Kalemci et al. 2006a, GX 339-4 during its decay in 2005 
Kalemci et al. 2006b, XTE J1720-318 CadoUe Bel et al. 
2004). Miyakawa et al. (2008) investigated the presence 
of cut-off from all bright hard state observations of GX 
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339-4 observed with HEXTE on RXTE, yet the statis- 
tics were not good enough to constrain the evohition 
of the cut-ofF parameters. Utihzing INTEGRAL ISGRI 
data in addition to the HEXTE data provided informa- 
tion about the evohition of the cut-off from the hard 
state to the hard intermediate state in the rising phase 
of the outburst of GX 339-4 (Motta et al 2009). How- 
ever, even in the brighter outburst rise, the relation be- 
tween high energy cut-off parameters and presence of 
jets is not weh estabhshed. For example, with the same 
data set but with a slightly different set of instruments, 
Caballero Garcia et al. (2007) and Joinet et al. (2008) 
reached conflicting results about the presence of cut-offs 
in the hard state spectrum of GRO J1655-40 during the 
outburst rise while the jet emission was present. 

Measuring distances of Galactic black hole sources pro- 
vides important information about the birthplaces of X- 
ray binaries, and, more importantly, allows us to con- 
strain the luminosities of these sources. Obtaining the X- 
ray luminosity of these sources would allow us to tie spec- 
tral states to physical processes in the accretion physics. 
It was shown by Maccarone (2003) that transition lu- 
minosities occur at similar values. This can be related 
to the changes from an accretion regime dominated by 
the disk to one dominated by a corona. Furthermore, 
this could also help in determining if there is a critical 
luminosity at which jets can form, and if jet dominated 
states exist at low luminosities (e. g. Russell et al. 2010), 
at what luminosity this takes place. Moreover, luminosi- 
ties are important for studying the fundamental plane 
that connects the radio - X-ray correlation Corbel et al. 
(2000, 2003) from stellar mass black holes to the su- 
permassive black holes (Merloni et al. 2003; Falcke et al. 
2004; Kording et al. 2006) as well as the variability stud- 
ies that connect these black hole sources (Kording et al. 
2007). 

1.1. XTE J1752-223 

XTE J1752-223 was discovered in the Galac- 
tic bulge region with RXTE on 23 October 2009 
(Markwardt et al. 2009b). It was suggested to be 
a black hole candidate by Markwardt et al. (2009a). 
Strong, relativistic iron emission lines are detected by 
Suzaku and XMM-Newton (Reis et al. 2011). The source 
showed typical outburst evolution of GBHTs, and it 
was intensely monitored using several satellites: RXTE 
(Muiioz-Darias et al. 2010; Shaposhnikov 2010), MAXI 
(Nakahira et al. 2010) and Swift (Curran et al. 2011). 
A distance of 3.5±0.4 kpc and black hole mass of 8- 
IIM0 were determined by Shaposhnikov et al. (2010) 
using QPO frequency saturation and comparison with 
other sources. Ratti et al. (2012) recently discussed the 
X-ray - radio correlation in this source. They also sug- 
gested that the source is in the Galactic bulge or closer 
to us, i.e. < 8 kpc based on comparison of transition 
luminosities of black hole transients (Maccarone 2003). 
The source's core location was accurately determined us- 
ing the astrometric optical observations and the VLBI 
radio imaging (Miller- Jones et al. 2011) in addition to 
the detection of decelerating jets and receding cjceta 
(Yanget al. 2011). 

Russell et al. (2012) discussed a late jet re-brightening 
in the HS via multiwavclength observations during its 
outburst decay towards quiescence (this is called flare 3 
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Figure 1. Hardness-intensity diagram of XTE J1752-223. Only 
RXTE PGA PCU2 data were used. Y-axis is count rate in 2.8 - 25 
keV. Hardness is defined as the ratio of (5.2 - 9.8 keV) to (2.8 - 5.2 
keV) band count rates. The filled circles are the observations used 
in this work. The vertical arrow show the hardness/intensity at 
the timing transition. The dotted line show the hardness/intensity 
during the INTEGRAL observation. The dashd line show the hard- 
ness/intensity at the VLBI detecton of the core. The dot-dash line 
show the hardness/intensity at the near infrared peak when the 
ATCA radio spectrum is consistent with optically thick emission. 

in our paper, see Appendix). They suggested that the 
brightening in the optical and X-rays probably has a 
common origin, but it is not clear if the origin is syn- 
chrotron radiation from a compact jet. They also sug- 
gested that the jet power changes at the peak of the flare 
and that there is evidence that the jet break between 
optically thin and thick synchrotron emission shifts to 
higher frequencies during this flare. 

In this paper, we present an in-depth multiwavc- 
length analysis of XTE J 1752-223 during its outburst 
decay, covered well with RXTE and Swift in X-rays, 
SMARTS in 0/NIR, and Australia Telescope Compact 
Array (ATCA) in radio. The spectral and temporal 
analyses for the source were examined for the whole 
outburst decay by simultaneous monitoring with RXTE 
and Swift during MJD 55,240-55,390. Fig. 1 shows the 
hardness-intensity diagram of the source during the en- 
tire outburst obtained using the RXTE data. We re- 
port on the observations for which simultaneous Swift 
data are available which are shown with flUed circles in 
Fig. 1. The X-ray behavior of the source is compared 
to the evolution in NIR and optical as well as the ra- 
dio. The decaying part of the outburst was covered 
previously with multiwavclength studies by Russell et al. 
(2012) and Brocksopp et al. (2013). However, our work 
enhances previous studies in several ways. We show for 
the first time the SMARTS light curve in the optical and 
the near infrared, which allows us to investigate the rela- 
tion between the radio emission. X-ray emission and opti- 
cal/NIR emission in much more detail than before. In ad- 
dition, we use RXTE and Swift observations together for 
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Figure 2. The evolution of spectral and temporal properties of 
XTE J1752-223. The photon index was obtained from combined X- 
ray spectra of XRT and PCA. The rms and flux after MJD 55,303 
were corrected due to the ridge emission in PCA. The fast transi- 
tion determined by timing analysis is indicated by a dashed line at 
MJD 55,282. The dotted Une shows the time of the INTEGRAL 
observation. The two panels on the bottom are of the / and H 
band NIR light curves, respectively. The observations shown in 
diamonds in the top panel correspond to those for which we show 
the power spectra in Fig. 4 

X-ray fits, and model the effects of Galactic ridge emis- 
sion to obtain reliable spectral and temporal parameters 
at low flux levels. The X-ray timing properties of the 
source and their relation to the multiwavclcgth emission 
arc investigated in terms of both evolution of rms am- 
plitudes and the shape of the power spectra. Moreover, 
an additional 3-day (MJD 55,305-55,307) broadband X- 
ray spectrum {Swift, RXTE, and INTEGRAL), covering 
0.6 - 200 keV, was also obtained in order to constrain 
the high energy behavior of the source (the time of the 
INTEGRAL obsrvation is shown in Figs. 1, 4, and 3). 
Finally, we analyze the XMM-Newton RGS spectra to de- 
termine the Hydrogen column density (Nh) of the source 
and use this information to determine the distance. We 
compare the result with the previously reported distance 
estimates (Shaposhnikov et al. 2010; Ratti ct al. 2012). 

2. DATA REDUCTION 
2.1. RXTE 



The X-ray evolution of XTE J1752-223 in 2010 out- 
burst is shown in the top three panels of Fig. 2: unab- 
sorbcd flux in 3-25 keV, rms amplitude of variability, and 
power-law index, respectively. The source was monitored 
throughout the outburst decay, including the fast tran- 
sition from SS to HS at MJD 55,282 until it went below 
detection after MJD 55,380. 

The RXTE PCA data were reduced with the scripts 
developed at UC San Diego and University of Tiibingen 
using HEASOFT v6.7. In the extraction of the en- 
ergy spectra, the photons from all available PCUs were 
considered. The background spectra were created from 
"bright" or "faint" models on the basis of the net count 
rate being greater or less than 70 c/s/PCU, respectively. 
We added 0.5% systematic error to all PCA spectra as 
suggested by the RXTE team. 

We used Tiibingen Timing Tools in IDL to compute 
the power density spectra (PDS) of all observations using 
PCA light curves in 3-25 keV band. The dead time ef- 
fects were removed according to Zhang et al. (1995) with 
a dead-time of 10 ^s per event, and the PDS is normalized 
according to Miyamoto & Kitamoto (1989). Broad and 
narrow Lorentzians (quality factor greater than 2 are de- 
noted as QPOs) are used for fitting (Kalemci ct al. 2005; 
Pottschmidt 2002). The rms amplitudes are calculated 
by integrating rms-normalized PDS from Hz to infin- 
ity. The peak frequencies are calculated as described in 
Belloni et al. (2002). 

2.2. INTEGRAL 

We observed the source using INTEGRAL (the INTer- 
national Gamma- Ray Astrophysics Laboratory), obser- 
vation ID "07400260001" for dates between MJD 55,304 
and 55,306 (Rev. 0917, 0918). The IBIS and SPI data 
were reduced and analyzed with the help of the ISDC 
Off-line Scientific Analysis (OSA) software package 9.0, 
and only the ISGRI spectra were used for further analy- 
sis in order to get better statistics at the relevant energy 
range. The INTEGRAL data were divided into three 
parts with approximately equal observing time to check 
whether the spectrum evolves at high energies within 2 
days. Since the fit parameters did not vary within the 
errors, we decided to merge two revolutions in order to 
increase statistics, i.e., the total integration time became 
177.4 ks. The observation time is indicated by a dotted 
line in Fig. 2. Note that neither spectral nor temporal 
properties show significant evolution around the time of 
the INTEGRAL observation. 

2.3. Swift 

The Swift XRT data were processed with the standard 
procedures (xrtpipcline vO.12.3). In the observations we 
analyzed. Windowed Timing mode data were available 
before MJD 55,367 and Photon Counting mode data 
were used afterwards. The energy spectra of the source 
were extracted from the photons within a circle of radius 
35" centered at the source position. The background 
spectra were extracted from photons within regions be- 
tween 70" and 100" from XTE J1752-223 centroid. The 
selection of event grades was 0-2 for Windowed Tim- 
ing mode, and 0-12 for the Photon Counting mode. 
XRT response matrices swxwt0to2s6_20070901v012.rmf 
and swxpc0tol2s6_20070901v011.rmf were used for the 
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Table 1 

Log of X-ray spectral and temporal parameters 



MJD 


XTE Obsid 


Swift Obsld^' 


Flux'' 


rms (%) 


r 


T,„{keV) 


i/i (Hz)" note 


55242 


95360-01-04-04 


31532014 


26.328 






9 974-0.04 


0.57ini 






55252 


95360-01-05-02 


31532021 


14.236 






2 04+0-08 


0.53«;»1 






55257 


95360-01-06-04 


31640001 


12.397 






9 n5+o-08 


0.52+0;gl 






55266 


95360-01-08-00 


31640003 


7.7581 






9 10+0.04 
^•^■'-0.04 


0.49l°ri 






55270 


95360-01-08-03 


31640005 


6.1194 






2.091H1 


0.48«;«1 






55277 


95360-01-09-03 


31640006 


4.6238 






2.ioJ:H^ 


0A8t°°l 






55580 


95360-01-09-05 


31640007 


4.2723 






1 qs+0.04 


4Q+0.01 

'J-^'^-o.oi 






55284 


95360-01-10-02 


31640008 


10.445 


16.36±1.24 


9 n4+0.03 


4t;+0-0i 


3.15±0.68 


55285 


95360-01-10-04 


31640009 


14.444 


18.72±1.15 


9 14+0.03 


(,41+0.01 
"^•^i-o.oi 


1.83±0.58 


55289 


95360-01-11-00 


31640010 


14.334 


20.09±0.68 


2 14+0-01 
^•^^-0.02 


42+0-01 


1.67±0.24 


55293 


95360-01-11-04 


31640011 


14.924 


25.00±4.15 


-, 0-1+0.01 
l-*l-0.02 


29+0''2 
^■^^-0.01 


0.55±0.22 


55297 


95360-01-12-01 


31640012 


14.129 


25.46±4.68 


1 64+0-01 


31+0-04 

U.Oi_Q Q4 


0.35±0.12 


55300 


95360-01-12-03 


31688001 


12.157 


26.17±4.17 


1 fifi+0.01 
l-"O-0.02 


r. 9t-+0.03 
"^■^^-0.02 


0.36±0.13 


55302 


95360-01-12-04 


31688003 


10.993 


26.63±7.61 


1 71+0.02 
^•'1-0.02 


n 2fi+o-03 


0.41±0.10 


55304 


95702-01-01-01 


31688006 


10.462 


28.40±4.72 


1 c:q+0.02 


29+°"^ 
^^■^^-0.03 


0.35±0.10 ISGRI 


55306 


95702-01-01-03 


31688008 


9.7074 


25.15±2.51 


1 60+0-03 
1-0U_Q 03 


9C-+O.IO 


0.41±0.17 ISGRI 


55307 


95702-01-01-04 


31688009 


9.8888 


27.30±6.00 


1 c:q+0.02 
l-^y_0.02 


n. 9C+O.O2 

'J-^°-0.02 


0.29±0.11 


55310 


95702-01-02-00 


31640013 


8.2724 


26.44±4.30 


1 fj4+0.03 
1-O4_0.03 


31+0-08 

•J-^i-o.os 


0.33±0.13 


55311 


95702-01-02-01 


31640014 


8.2815 


27.41±2.99 


-1 48+0.01 


7+0.04 
'^•^'-0.03 


0.41±0.13 


55313 


95702-01-02-03 


31688012 


7.3498 


25.40±4.51 


-T f!7+0.04 


n 9c;+0.07 
O.^t)_0.04 


0.43±0.26 


55313.5 


95702-01-02-04 


31688013 


7.3846 


24.41±2.35 


-1 CI +0.02 
1-^1_0.03 


r. 00+0. 16 
'^■■'■^-0.06 


0.32±0.09 


55320.5 


95702-01-03-04 


31688014 


5.6153 


27.69±3.95 


1 f!o+0.03 


No disk 


0.36±0.12 


55328 


95702-01-04-04 


31688016 


3.1901 


26.65±2.74 


1 trq+0-05 


needed 


0.29±0.11 a single 


55333 


95702-01-05-03 


31688017 


1.6355 


36.62±4.26 


1 fj4+0.05 
1-O4_0.06 




0.17±0.06 Lorentzian 


55336.5 


95702-01-05-06 


31688018 


1.3928 


35.26±4.47 


1 fj4+0.06 




0.15±0.05 enough 


55340.5 


95702-01-06-02 


31688019 


1.2990 


38.94±4.83 


1 f-o+0.07 




0.15±0.04 


55343 


95702-01-06-03 


31688020 


1.1554 


40.67±6.84 


1 cc + O.ll 

i-t'O-o.io 




0.22±0.10 


55348.5 


95702-01-07-05 


31688021 


2.0856 


30.56±7.86 


1 f-1+0.08 




<1.41 


55352.5 


95702-01-08-00 


31688022 


2.7439 


35.94±4.28 


1 cs+0.03 




0.25±0.05 


55356 


95702-01-08-02 


31688023 


3.0706 


31.65±2.90 


1 cc+O.OS 

1.0J_Q Q4 




0.19±0.07 


55368 


95702-01-10-01 


31688024 


2.2719 


35.59±9.11 


1 co+0.08 
l-0o_0.08 




0.20±0.12 


55371 


95702-01-10-02 


31688025 


1.9037 


31.11±6.84 


1 cn+o.oT 

1-0U_Q Qg 




O.lSitO.lO 


55377 


95702-01-11-01 


31688026 


1.1313 


27.21±5.00 


1 70+0-10 

i-'i'-o.io 




0.09±0.04 


55380 


95702-01-12-00 


31688027 


0.2859 




1 44+0-21 
1-^^-0.21 




4.18±1.74 



^ Spectral parameters are obtained from the joint RXTE PGA and Swift XRT data fits 

1* The unabsorbed flux values between 3-25 keV are in units of 10~10 ergs/cm^ /s. 

'^ The timing properties of observations. Those shown in bold are discussed further in text and in Fig 

for the peak frequency of the Lorentzian that peaks at lower frequencies. 
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Windowed Timing and Photon Counting modes, respec- 
tively. We also generated auxiliary response files with 
the HEASOFT tool xrtmkarf. In making the ARF files, 
we used an exposure map produced with the tool xrtex- 
pomap. There was also some pile-up before MJD 55,272 
when the count rates were over 100 c/s, and for these 
observations we excluded the peak region of radius 5" . 

2.4. XMM-Newton 

The XMM-Newton observation of the source was per- 
formed on 6 April 2010, MJD 55,292. The total ex- 
posure time of the observation was 41.8 ks. Details of 
this observation and results from the EPIC-pn detector 
has been discussed in Reis et al. (2011). We here con- 
centrate on the RGS (Reflection Grating Spectrometer, 
den Herder et al. 2001) data, which provide high resolu- 
tion soft X-ray spectra in the range of 5 to 35 A. Our 
aim is to model the absorption edges in the soft X-rays 



to determine the hydrogen column density along the line 
of sight towards XTE J 1752 -223. 

We extracted RGS spectra using the rgsproc within 
SAS vll.0.0 with the latest calibration files available as 
of Dec. 2011. We grouped the RGS spectra with spec- 
group tool so that each spectral channel will have at least 
100 counts and the instrumental resolution will not be 
oversampled more than a factor of three. 

Given the X-ray flux of the source at the time of the 
observation, we also checked the RGS observation against 
any affect of photon pile-up. For this purpose we used 
the fluxed RGS X-ray spectra, as given by the Browsing 
Interface for RGS data (BIRD), and checked the flux on 
each CCD against the limits given in the XMM-Newton 
User's Handbook' . This comparison showed us that even 
though there was a photon pile-up it was smaller than 
2%. This is also thanks to the fact that the observation 

http://xmni.esac .esa. int/external/xnini_user_support/docimientation/uhl 
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was performed with double-node readout for RGSl and 
RGS2. Such a smah amount of saturation would have 
only minimal, if any, effect on the depths of individual 
absorption edges. 

2.5. Optical & NIR observations 

In addition to the X-ray observations, we analyzed the 
0/NIR SMARTS data (/ and H band) obtained with the 
CTIO 1.3m telescope, during MJD 55,240-55,390 cover- 
ing 150 days of the outburst decay. With the help of 
IRAF V2.14, we conducted photometry using the PSF 
(Point Spread Function) fitting due to the crowded region 
towards the source (see Russell et al. 2012 for the high 
resolution NIR image). Then, we selected five compari- 
son stars in the FOV as non- variables and went through 
the standard routines to get the actual magnitudes (see 
Buxton ct al. 2012 for standard procedures). The mag- 
nitudes we obtained could be affected by a nearby source, 
especially in the H band due to the poor observational 
conditions and the resolution of the telescope. The re- 
sults are shown in Figs. 2 and 3. 




5240 5260 5280 5300 5320 5340 5360 5380 
Date (MJD-50,000 days) 



Figure 3. The evolution of SMARTS O/NIR magnitudes and 
radio fluxes during the decay. The solid lines indicate exponentially 
decaying emission possibly from the outer disk, which are shown 
to demonstrate the presence of the three flares mentioned in the 
text (see Appendix). The radio data are added with diamonds 
and triangles (see legend). The errors on radio before MJD 55,320 
are smaller than the symbol size. The arrow represents the upper 
limit of the VLBI core radio emission at MJD 55,277 (sec text). 
The dates of the fast transition and the ISGRI observation are 
indicated by dashed lines on MJD 55,282 and 55,306, respectively. 
A colored version of this figure is available in the online version. 



2.6. Radio 

The ATCA provided the evolution of the source 
throughout the decay in radio in two bands, 5.5 and 
9 GHz. The peak fluxes in two bands were obtained 
with the help of the radio interferometry data reduc- 
tion package, Miriad (Sault et al. 1995). Details of 
the analysis and an extended discussion of the behav- 
ior of ths source in radio for the entire outburst can 
be found in Brocksopp et al. (2013). In addition, very 
long baseline interferometry (VLBI) measurements were 
also conducted for a few observations during the de- 
cay with the European VLBI Network (EVN) and the 



Very Long Baseline Array. The high resolution of 
VLBI made it possible to resolve a compact jet core on 
MJD 55,311 and 55,315 (Miller- Jones et al. 2011) (X-ray 
hardness/intensity at the time of VLBI detections are 
shown in Fig 1). We used the radio flux values reported 
by Yang et al. (2011) for the core (denoted as component 
D in their paper). In Fig. 3, we also show the 3 a upper 
limit of Q.llmJy for the core on MJD 55,277 based on 
the EVN image sensitivity (1 sigma = 0.03 mJy/beam 
in Table 1 of Yang et al. 2011) and -80% (typical value) 
phase-referencing coherence loss (Jun Yang, personal 
communication) . 

3. ANALYSIS AND RESULTS 

All i?XTE and 5'wift spectra were fitted with a phe- 
nomenological model of disk blackbody (diskbb) plus 
power law with absorption. A smeared edge model 
(snicdge, Ebisawa et al. 1994) was also added to the 
RXTE spectrum to improve the x^ similar to previous 
work by this group (Kalemci et al. 2004, 2005, 2006a). 
For photoelectric absorption, we used Tuebingen-Boulder 
ISM absorption (TBabs) with the abundance values of 
Wilms et al. (2000) and the cross sections of Verner et al. 
(1996). This choice is discussed in more detail in Section 
3.4. In addition, due to the proximity of the source to 
the Galactic plane (6 = -1-2°. 11), the Galactic ridge con- 
tribution is significant for RXTE when the source is at 
low luminosity levels. To take into account the Galactic 
ridge emission, an additional power-law was added to the 
spectral fit after MJD 55,310, fixing the index to 2.1 and 
normalization to 1.209 x 10~^ based on values given in 
Revnivtsev (2003). bfThe ridge emission count rates ob- 
tained with this model is consistent with the background 
rates for the PGA Galactic bulge scans^.The rms am- 
plitudes were corrected for the ridge contribution after 
MJD 55,303 as well as for background emission following 
Berger & van der Klis (1994). 

3.1. Multiwavelength evolution 

Fig. 2 shows the evolution of the X-ray spectral and 
temporal properties as well as the evolution in O/NIR 
magnitudes. The vertical line on MJD 55,282 indicates 
the transition in timing due to the fast change in the 
rms. An increase in the power-law flux accompanies the 
increase in the rms amplitude of variability. The X-ray 
flux continued to increase for a couple of days after the 
transition, whilst the photon index remained the same 
until the peak in X-ray flux. As the X-ray flux started 
to decay, the photon index abruptly decreased to '^1.7. 
The spectrum continued to harden and the index leveled 
off around 1.6. A secondary peak was also observed in 
X-rays after MJD 55,340. 

The SMARTS optical and NIR observations showed in- 
teresting behavior as shown in Figs. 2 and 3. Some flares 
are present in the light curves at around MJD 55,280 
(flare 1), MJD 55,320 (flare 2) followed by a large increase 
in flux starting around MJD 55,335 (flare 3) coinciding 
with the increase in the X-ray flux (RusseU et al. 2012). 
The small increase during flare 1 in the H band coin- 
cides with the timing transition which is also reported 
in Buxton et al. (2010). We note that after analyzing 

* http : //asd . gsf c . nasa. gov/Craig . Markwardt/galscaii/html/XTE_ J1752-22. 
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the data in the fohowing days, it is found that the in- 
crease is not significantly above the continuum as shown 
in Fig. 3. Flare 2 starting at MJD 55,320 in the / band 
coincides with the end of hardening in the X-ray spectra. 
See Appendix for the discussion of significance of fiarcs. 
Fig. 3 also shows the evolution of radio fiux during out- 
burst decay at different wavelengths from different obser- 
vations. The green and red diamonds are from ATCA at 
5.5 GHz and 9 GHz, respectively. We also showed the 
H and / band magnitudes to associate changes in radio 
to the changes in optical/NIR. Until MJD 55,320, the 
radio spectrum is optically thin. At MJD 55,311 and 
MJD 55,315, the radio core is detected with the VLBI, 
while the ATCA radio spectrum was still optically thin. 
This may indicate that there is contamination in the ra- 
dio data, and the ATCA flux includes emission from a 
compact jet and some other interaction. The detection 
of the core corresponds to a little flare that can be seen 
in the / band. Around MJD 55,335 the ATCA spectrum 
gets harder, and this is where the fiux in /, H band, and 
X-rays increase. 

3.2. Timing evolution 

As distinct multiwavelength changes occur during the 
brightening after MJD 55,340, we decided to check 
whether there are any changes in the PSD at the same 
time. The evolution of PSDs is shown in Fig. 4. We 
started from day MJD 55,285, during the first X-ray 
peak, and went down in fiux to compare observations 
before and during the peak with similar fiuxes. At 
MJD 55,285 two Lorentzians are required to fit the data, 
and the first Lorentzian peaks at '^1.83 Hz. The rms 
amplitude of variability is --18.72%. At MJD 55,297 the 
spectrum is now hard (r^l.7), and the rms amplitude is 
higher (^25.46%). Still, two Lorentzians are required to 
fit the data, but the peak frequency of the first Lorentzian 
is shifted down to ^0.35 Hz. For the rest shown in Fig. 4, 
a single Lorentzian is enough to fit the data. We note 
that the ridge correction significantly affects the rms am- 
plitudes for flux levels less than 3 x 10"^" ergs cm~^ 
s~^, which is reflected as large errors. At MJD 55,328 
and MJD 55,356, the X-ray fluxes are almost the same 
with similar spectral parameters. The PSDs are also 
similar, and the rms and peak frequencies are consis- 
tent within errors. Likewise, the PSDs at MJD 55,343 
and MJD 55,371 are also at similar flux levels, and show 
no signiflcant change in spectral and timing properties 
(see Table 1 for parameters). Even though the statistics 
are not good enough to demonstrate a significant differ- 
ence (or lack thereof), the timing properties appear not 
to change before and during the flare, and similar X-ray 
spectral properties provide similar PSDs. 

3.3. Broadband X-ray spectrum 

To test whether spectral breaks disappear with the 
formation of jets, we conducted observations with the 
INTEGRAL observatory. The ISGRI spectrum is com- 
bined with contemporaneous RXTE PCA data. We 
started with our standard model, consisting of absorp- 
tion, diskbb, smedge and power-law. We saw strong 
residuals indicating a high energy cut-off component in 
the fit (see Fig. 5a). Adding a high energy cut-off (high- 
ecut) significantly improved the fit (F-tcst results in a 
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Figure 4. Several selected PSDs from the intermediate state to 
low flux hard state. The vertical axes are frequency X powrer and 
the dates (MJD-50,000) are indicated on the left of each panel. 
The peak of the Lorentzian gradually decreases. The bottom three 
panels correspond to times before, during, and after the largest 
NIR flare, and there are no clear differences between these power 
spectra. 

chance probability of 10~^", see Fig. 5b and c). The fit 
parameters can be found in Table 2. 

Moreover, we also tried physical models and tried to 
establish whether an iron line is present in the spec- 
trum. We fit the combined ISGRI, PCA, XRT spectrum 
first with the thermal Comptonization code, CompPS 
(Poutanen & Svensson 1996, see Fig. 5d) and obtained a 
good fit with a reduced x^ of 0.97. The statistical qual- 
ity of the XRT spectrum was not good at the iron line 
region. We added diskline (Fabian et al. 1989) to model 
the iron line emission at '^6.4 keV. The improvement in 
the fit was not significant. Our dataset does not allow 
us to constrain the iron line parameters of this source at 
the time of the INTEGRAL observation. 

3.4. Distance to the source 

Near-IR observations of red clump stars have been 
used to map the extinction as a function of distance 
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Figure 5. Spectral fits to the Swift XRT (shown with + sign, colored black in the online version), RXTE PCA (shown with diamonds, 
colored purple in the online version) and INTEGRAL ISGRI data (shown with triangles, colored red in the online version), a) Count 
spectrum and fit with a power-law and disk blackbody. b) Count spectrum and fit with a power-law and a high energy cut-off and disk 
blackbody. c) Unfolded spectrum with high energy cut-off. d) Unfolded spectrum with Comptonization model compps. A colored version 
of this figure is available in the online version. 



Table 2 

The best fit parameters for simultaneous XRT, PCA, and ISGRI spectra 



Models 



Nh 



kTf 



bb, 



''II^T," 



kT^{keVY 



tikeVr Ef.uikeVY xVidof) 



CompPS+diskline 
power*highecut 0.67' 



U.DU_Q Q2 

-1-0.04 

0.04 



7-f0.03 
-0.02 



0.27"^ 



1.24 



-1-0.24 



0.53; 



105.2 



-1-14.8 



1.67; 



-0.02 
-0.02 



25.4"* 



236.5 



312(323) 



+''^* 342(323) 



"• The inner radius temperature of soft photons in multicolor disk. 
^ The vertical optical depth of the corona. 
'^ The reflection factor. 
"* The electron temperature in the corona. 
^ The cut off energy. 
The folding energy. 

along a given line of sight in the Galaxy, as reported 
in a number of studies by Paczynski & Stanek (1998); 
Lopcz-Corrcdoira et al. (2002); Cabrcra-Lavers et al. 
(2005); Nishiyama et aL (2006). As they have a very nar- 
row luminosity function (especially in the near-infrared) , 
and they can be easily isolated from a color magni- 
tude diagram, they have been largely used as stan- 
dard candles in describing the geometry in the inner 
Galaxy (Cabrcra-Lavers et al. 2007, 2008). For all the 
above, this population becomes a very powerful tool both 



for tracing and modeling the interstellar extinction in 
the more obscured areas of the Milky Way (see, e.g., 
Drimmel et al. (2003) or Marshall et al. (2006)). 

These stars have also been used to estimate the dis- 
tances of X-ray sources by comparing the estimates 
for the increase in extinction with the distance derived 
from the NIR data with the intrinsic extinction of the 
source derived from their X-ray data, with very suc- 
cessful results (see, e.g., Durant & van Kerkwijk (2006); 
Castro-Tirado et al. (2008); Giiver et al. (2010)). Here, 
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Table 3 

Resulting Hydrogen column density values obtained from fits to individual absorption edges in 
the RGS data. Values are in units of 10^^ em^^ 



O 



Fe 



Ne 



Mg 



Average'' 



ISM"" 
Solar == 



1.04±0.05 
1.04±0.05 



0.88±0.04 
0.73±0.04 



1.17±0.12 
1.17±0.12 



1.08±0.18 
0.99±0.14 



0.96±0.03 
0.88±0.03 



^ Error wreighted averages of the values found from individual edges. 

*> As given by Wilms et al. (2000). 

'^ As presented by Asplund et al. (2009). 




i%#^ 



^ 



^ 



^i^^^^^\^^^VT^^'^mi^-^ '^ - ^ 



#i 



Wavelength ( A ) 




^^++***f¥+*=F%+ 




^ + 



\ii^^i 



_° L^++^44^^v+^r^^H#^++vv ^ 



Wavelenglh ( A ) 




Wavelength ( A ) 



Wavelength ( A ) 



Figure 6. RGSl (blaek) and RGS2 (red) data at the Mg (upper left), Ne (upper right), Fe (lower left), and O (lower right panel) absorption 
edge regions arc shown in the upper panels of each plot. Lower panels show the residuals from the best fit model. Sharp features are due 
to CCD gaps in the RGS detector. A colored version of this figure is available in the online version. 



we follow a similar method to map extinction vs. dis- 
tance in the direction of XTE J1752-223 by using pho- 
tometric data from the 2MASS survey (Skrutskie et al. 
2006). 

To find the intrinsic extinction in NIR, we first mea- 
sured the interstellar X-ray absorption of the source in 
a way that is independent of assumptions about the in- 
trinsic X-ray spectral properties by modeling individual 
absorption edges of the elements O, Ne, Fe, and Mg. For 
this purpose we only use the X-ray grating observation 
of the source obtained with the RGS on board XMM. 

To measure the Hydrogen column density towards 
XTE J1752-223, we followed a method similar to 
Durant & van Kerkwijk (2006) and Giiver et al. (2010). 
We selected only a small wavelength range for each 
edge; 8.5-10.5 A region for Mg, 13.5-15.0 A for Ne, 
16.0-18.0 A for Fe, and 20.0-25.0 A for the O edge 
(note that in this wavelength range, RGS2 has no sen- 
sitivity because of a CCD failure early in the mission. 
Hence, we did not include data from RGS2 for the O 



edge). We assumed that the continuum in these small 
intervals can be modeled with a power-law function and 
modeled each edge using the tbnevP model (Wilms et al. 
2011, in prep). For spectral modeling we used XSPEC 
12.7.0 (Arnaud 1996). We used the elemental photoelec- 
tric absorption cross-sections as given by Verner et al. 
(1996). Finally, we performed the fits both assuming ISM 
and solar abundances as given by Wilms et al. (2000) 
and Asplund et al. (2009), respectively. We present the 
resulting Hydrogen column density values in Table 3 for 
each element and show our fits to the data in Fig 6. 

In the O edge region (at 23.05 A) we noticed an absorp- 
tion line like feature in the RGSl spectrum that is not 
modeled in the tbnew model. We modelled this feature 
with a Gaussian to obtain a better fit for the absorp- 
tion edge. We note that adding this component does not 
change the measured O column density although adding 
it improves our statistics. The position is consistent with 

^ http://pulsar.sternwarte.uni-erlangen.de/wilms/researeh/tba- 
bs/ 
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Figure 7. Left: Noar-IR color-magnitude diagram for the field centered around XTE J1752-223. Red crosses show points of maximum 
density of red clump stars for each individual magnitude bin and their uncertainties. Right: Evolution of near-IR extinction in the direction 
of XTE J1752-223 as derived from the 2MASS archive. A colored version of this figure is available in the online version. 



S XIV resonance line. Similarly, we noticed a statistically 
significant absorption line like feature in the Ne edge re- 
gion both in the RGSl and RGS2 data at 14.6 Aas well. 
Again we modelled this feature with a Gaussian to obtain 
a better fit for the continuum. This position is consistent 
with Fe XVIII resonance line and the line can be seen in 
Fig. 6. 

Following the method outlined in Giiver et al. (2010) 
and the relations presented by Giiver & Ozel (2009) 
and Cardelli et al. (1989), the Hydrogen column den- 
sity we measure assuming an ISM abundance results 
in an optical extinction of Av=4.34±0.22 mag or near- 
infrared extinction of Akj.= 0.49±0.09 mag. Similarly 
assuming solar abundances, we get an optical extinc- 
tion of Av=3.98±0.21 or near-infrared extinction of 
Ak,=0.46±0.08 mag. 

In Fig. 7 (right), we show the derived evolution of 
the extinction as we move toward the inner Galaxy in 
the direction of XTE J1752-223 {1^6° A2 , 6=+2°.ll), 
while the region formed by the red clump stars can be 
seen in the color- magnitude diagram (see Fig. 7, left). 
It can be seen that the extinction increases up to 4.5-5 
kpc due to the interstellar dust in the Galactic plane in 
this direction, and remains nearly constant for distances 
larger than these. We observe this flattening since we 
have reached the bulge component at this distance along 
the line of sight, making it harder to accurately trace the 
extinction for larger distances. The red clump sources 
for the Bulge at a Galactic latitude of 2° above the plane 
completely dominates the star counts at this distance, 
and the (few) red clump sources for the Galactic disc 
are completely diluted. Hence, extinction at greater dis- 
tances is uncertain. Therefore, this method of using 
red clump stars, unfortunately, has limited capability 
in terms of measuring the distance to XTE J1752-223. 
However, based on the extinction values we obtained for 
the source from X-ray spectroscopy and assuming that 
at least a significant part of the X-ray absorption is not 
due to intrinsic absorption, XTE J 1752-223 cannot be 
closer than 5 kpc. 

4. DISCUSSION 



4.1. State transitions and jets 

The results of X-ray monitoring and the 0/NIR light 
curves are shown in Fig. 2. and the radio evolution is 
included in Fig. 3 (see also Brocksopp et al. 2013). Flare 
1 in the H band around MJD 55,280 is not significant, 
yet, it corresponds to a change in X-ray timing properties 
(Buxton et al. 2010), indicating a transition to an inter- 
mediate state, as well as an increase in the ATCA radio 
flux between MJD 55,275 and MJD 55,290 denoted as 
"Peak 7" in Brocksopp et al. (2013). Close to the begin- 
ning of this flare, on MJD 55,277, the core is not detected 
with the VLB I with an upper limit of 0.11 mJy. Accord- 
ing to Brocksopp et al. (2013), the increase in the radio 
flux is probably from the interactions of jet components 
rather than compact jet turning on. 

On the other hand, flare 2, between MJD 55,310 and 
MJD 55,330, is coincident with both the VLBI detec- 
tion of the core on MJD 55,312 and a small increase 
of the ATCA flux. Since flare 2 also coincides with 
the end of hardening of X-ray spectra, it is more likely 
to be a result of compact jet forming (Dinger et al. 
2012; Kalemci et al. 2013), supporting the conclusion 
of Brocksopp et al. (2013). Unfortunately, the lack of 
VLBI coverage between MJD 55,277 and MJD 55,312 
prevents us from coming to strong conclusions on when 
exactly the compact core turned on. 

The ATCA spectral index increases towards 
MJD 55,350, and the spectrum becomes consistent with 
a flat spectrum on MJD 55,350 (hardness/intensity 
at the time are shown in Fig. 1). This spectral 
evolution coincides with the flare 3 seen in both the / 
and the H band starting at MJD 55,335. A compact jet 
may be dominating the 0/NIR and radio emission after 
this date (also see Russell et al. 2012). An increase in 
0/NIR flux along with a radio spectral evolution from 
an optically thin to optically thick emission is also seen 
for GX 339-4 (Corbel et al. 2012). 

If the emission from a compact jet dominates the 
X-ray emission during the second long flare between 
MJD 55,340 and 55,380 (Russefl et al. 2012), one could 
expect some change in the timing properties before and 
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after the peak. We did not observe a significant change 
in the PSDs, their rms amplitudes, and the peak frequen- 
cies when the compact jet is present. However, because 
the errors of the peak frequencies in these datasets are 
not neghgible, ^30%, timing studies for brighter sources 
might be more useful to discuss the effects of jets on 
timing. 

4.2. Spectral breaks in hard X-rays 

The presence of spectral breaks is well established 
during the transition to (or from) the hard state 
(Kalemci et al. 2005; Motta et al. 2009). There is some 
indication that the breaks disappear after the jet turns 
on based on spectral fitting of data from HEXTE 
(Kalemci et al. 2006b). However, there is no conclusive 
evidence that relates the jet formation to the evolution 
of spectral breaks as the short monitoring observations 
with RXTE often lack the statistics required to charac- 
terize this evolution. To obtain better statistics in hard 
X-rays, we triggered our INTEGRAL observation around 
the time that the spectral index is hardest. The obser- 
vation took place just before the core is detected wth 
the VLBI. The broadband combined spectrum clearly 
shows a break, with a folding energy at ^237 keV. This 
result is consistent with earlier work showing that a cut- 
off is present before a compact jet dominates the radio 
and 0/NIR emission. Unfortunately, the timing of the 
observation did not allow us to test the hypothesis that 
the cut-offs disappear after the strong compact jets dom- 
inates the multiwavelength emission. The hard X-ray 
spectrum of the source is consistent with thermal Comp- 
tonization. 

4.3. Distanee 

The method using red clump stars provides a lower 
limit to the distance to the source of > 5 kpc. 
Shaposhnikov et al. (2010) estimate the distance of the 
source as about 3.5 kpc. At this distance, the Galactic 
extinction is about 60% of what we infer from the ab- 
sorption edges observed with XMM RGS. Assuming that 
the intrinsic absorption in this source is not as high as 
2 magnitudes in the optical, this distance estimate is in- 
compatible with our results. Such a high intrinsic absorp- 
tion seems hard to understand given that no dips in the 
X-rays have been observed, which could have been inter- 
preted as neutral matter absorbing X-rays emitted from 
the boundary layer. It has been found by Miller et al. 
(2009) that the X-ray absorption as derived from pho- 
toelectric absorption edges remains constant as the lu- 
minosity and spectral states of X-ray binaries, including 
black hole systems, vary. This finding suggests that the 
ISM strongly dominates the measured neutral Hydrogen 
column density in the spectra of X-ray binaries. 

Also, possible uncertainties in the extinction law 
used for deriving the interstellar extinction (see, 
e.g., Nishiyama et al. 2009; Gonzalez- Fernandez 
2012, in press) cannot support such high differences 
from the estimate by Shaposhnikov et al. (2010), as 
the affects of such differences are within the range 
covered by the errors shown in Fig. 7. Hence, the 
assumption of XTE J 1752-223 being no closer than 
5 kpc is not produced by uncertainties in the NIR 
extinction measurements. 



We note that the distance of 3.5 kpc is also not 
compatible with the overall behavior of GBHTs in 
terms of their transition luminosities (Maccarone 2003; 
Ratti et al. 2012). As shown in Kalemci et al. (2013), 
the overall luminosity evolution of this source becomes 
compatible with other black hole transients if we assume 
a distance around 8 kpc, and a distance of 3.5 kpc would 
make the behavior of this source quite different from the 
other black hole binary systems. Future observations of 
this source during a new outburst can test our measure- 
ments and provide clues about the any variabilities in 
the amount of X-ray absorption, which would be a direct 
evidence for intrinsic absorption. 

Apart from the evidence from the transition luminosi- 
ties, the radio - X-ray correlation is also consistent with 
the behavior of other radio-quiet sources if the distance 
is 8 kpc (Brocksopp et al. 2013). The larger distance not 
only boosts the minimum power in the jet, but also in- 
creases the speed of the jet. As indicated by Yang et al. 
(2011), a larger distance means that the XTE J1752-223 
may be a superluminal source. 

5. SUMMARY AND CONCLUSIONS 

In this work, we analyzed data from the Galactic black 
hole binary XTE J 1752 -223 from the radio band ah the 
way to hard X-rays during its outburst decay in 2010. We 
investigated the evolution of X-ray spectral and tempo- 
ral properties of the source using RXTE and Swift and 
compared the results to the evolution of fluxes in the 
/ and H band SMARTS data and also to the spectral 
evolution of the ATCA radio data. We also studied the 
XMM-Newton RGS data to find the extinction towards 
the source. The important results from this analysis are 
summarized below: 

• We showed that the detection of the radio core with 
the VLBI corresponds to a small fiare in the / 
band flux (flare 2 in Fig. 3). This flare is followed 
by a larger flare in the / and H band flux (fiare 3 
in Fig. 3). 

• Flare 3 in the / and the H band coincides with the 
ATCA radio data becoming optically thick, relat- 
ing the compact jet to the 0/NIR changes. 

• The short term X-ray timing properties do not 
show a significant change during fiares in the 
0/NIR. 

• The broadband X-ray spectrum including 
INTEGRAL ISGRI indicates a high energy 
cut-off, which is also seen in other sources before 
strong compact jets are observed. The spectrum 
is consistent with thermal Comptonization. 

• We showed that the distance to the source cannot 
be less than 5 kpc. Based on the analysis of transi- 
tion luminosities and also the behavior of of radio 
- X-ray flux correlation, the distance is more likely 
to be close to 8 kpc. 
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Figure 8. Straight line + Gaussian fits to the flares discussed in the article. Top panel shows the / band light curve and the bottom panel 
shows the H band light curve. The dotted lines (colored purple in the online version) are linear fits to the data. The long dashed curves 
(colored red in the online version) are the fits with straight line+Gaussian for the flare 3. The dotted-dashed curves (colored orange in 
the online version) are fits with straight line + 2 Gaussians. The solid curves (colored black in the online version) are the overall fits with 
straight line + 3 Gaussians. The solid, straight line represents the linear fits to the continuum when all Gaussians are present. A colored 
version of this figure is available in the online version. 



APPENDIX 
SIGNIFICANCE OF FLARES 

We have discussed three flares shown in Fig. 3. The first flare was the subject of Buxton et al. (2010), and the 
third flare was discussed in detail in Russell et al. (2012). Close to the dates of the detection of the radio core, we 
realized that there is some excess emission especially in the / band. We then tested the significance of all flares in 
the / and the H band. To do that, we assumed that there is underlying emission that exponentially decays, which 
would correspond to a linear decrease in magnitudes (Buxton & Bailyn 2004; Kalcmci et al. 2005; Russell et al. 2010; 
Buxton et al. 2012; Dinger et al. 2012). We first fitted the entire datasets with straight lines (dotted lines in Fig. 8). 
For both bands the fits were very poor (the x^ values of all the fits are given in Table 4). Following Buxton & Bailyn 
(2004), we then fitted the excesses over the continuum with Gaussians. 

We started with the most obvious excess on '^MJD 55,370 (flare 3, long dashed lines in Fig. 8). This flare is required 
in both bands. Some residuals were present in the / band around MJD 55,320 and around MJD 55,290 in the H 
band. Adding Gaussians in the fit to these flares (dotted-dashed lines in Fig. 8) we found that flare 2 is statistically 
significant in the I band whereas flare 1 in the H band is not statistically significant. Finally, for completeness, we 
added Gaussians around MJD 55,290 in the / band and around MJD 55,320 in the H band and refit the datasets 
(solid curves in Fig. 8). These Gaussians are not formally required in the fit (chance probabilities obtained using the 
F — test are given in Table 4) . 
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Table 4 

Fit results to the O/NIR light curves 



Flare # 


xf'^ DOFi'' xi'^ Z)OF2<* 


pc 




Fits to /-band light curve. 




3 


163.40 62 80.72 59 


4.1410-9 


2 


80.72 59 58.77 56 


4.5610-* 


1 


58.77 56 50.44 53 
Fits to JZ-band light curve. 


0.043 


3 


63.72 55 40.41 52 


2.6210-5 


1 


40.41 52 35.91 49 


0.12 


2 


35.91 49 35.50 46 


0.91' 



Degrees of freedom before adding the new component. 
'^ X^ after adding the new component. 
'^ Degrees of freedom after adding the new component. 
" Chance probability from the F-distribution. 

To get a fit, the width of the Gaussian was limited to 3-5 days. 



